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^— ^ . Abstract 

, In this paper, the asymptotic formulas for Eulerian numbers, refined 

Eulerian numbers and the coefficients of descent polynomials are obtained 
directly from the spline interpretations of these numbers. Having related 
these numbers directly to B-splines [15], we can take advantage of many 
• powerful spline techniques to derive various results of these numbers. The 

asymptotic formulas for the Eulerian numbers Ad,k agree with the pre- 
viously known results which were given by L. Carlitz et al.(1972)[2] and 
S.Tanny (1973) [IB], but the convergence order is much better. We also 
give the asymptotic representations of refined Eulerian numbers which is 
in terms of the Hermite polynomials. 
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O ' 1 Introduction 

O ' 

t-H . 

Eulerian number, denoted here by Ad t k is the number of permutations in the 

•l-H , 

, symmetric group Sd that have exactly k—1 descents. They play important roles 

d ' in enumerative combinatorics, geometry [T^], statistical applications 2, l<f], and 

spline theory [TU1 115) . The asymptotic properties of Eulerian numbers was ex- 
amined by S. Tanny [TB] and L. Carlitz et. al. |2] with the help of central Limit 
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theorem of probability theory. They showed that the Eulerian numbers approx- 
imate to the Gaussian function with the convergence order 0(d~i). Neverthe- 
less, in this paper, we obtain a much better convergence order 0(d~i) with the 
help of the spline theory. 

The refined Eulerian number A-d,k,j is the number of permutations in the 
symmetric group Sd with k descents and ending with the element j. Brenti and 
Welker [5] gave a precise description of the h-vector of the barycentric subdi- 
vision of a simplicial complex in terms of the h-vector of the original complex 
and the refined Eulerian numbers. Later, Martina Kubitzke and Eran Nevo 
|13j derived new inequalities on A^fej by the algebraic method. A result due 
to Ehrcnborg, Readdy and Steingn'msson [T7J gave a geometry interpolation of 
A-d,k.j as the mixed volumes of two adjacent slices from the unit cube. In |15j 
Renhong,W., Yan, X. and Zhiqiang, X. gave a spline interpolation of the Ad.k.j 
by the geometry interpolation and the spline theory. Further more, they gave 
the explicit and recurrence representation of Ad^kj- In [Hj . the experimen- 
tal evidence suggested that the refined Eulerian number Ad,k,j is not exactly 
normal approximation as d grows large. In this paper, we give the asymptotic 
representations of refined Eulerian numbers in terms of the Hermite polynomi- 
als. 

Recently, Steingn'msson [7J generalized the definition of descents and ex- 
ccdanccs to the elements (called indexed permutations, colored permutations or 
r-signed permutations see [TTJ 13 03 E] ) of groups S% ■— Z n lSd, where I is wreath 
product with respect to the usual action of Sd by permutations of [d] . The num- 
ber of indexed permutations in S"2 with k descents is denoted by D(d,n, k). In 
[7J, using the work of Brenti [5], Steingrimsson [7J showed that the numbers 
D(d,n,k) are unimodal. Moreover, In|15j. using the spline theory, the authors 
gave the explicit expression of D(d, n, k) and showed the numbers D(d, n, k) are 
log-concave. In this paper, we show that D(d, n, k) is approximately normal in 
terms of the spline interpretation of D(d, n, k). 
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The paper is organized as follows. In Section 2, after recalling some neces- 
sary definitions and notations, we show the connection between the B-splines 
and these combinatorial numbers. In Section 3, we give the main results of 
this paper. The results show the development of various asymptotic properties 
of Eulerian numbers, refined Eulerian numbers and the coefficients of descent 
polynomials, which are obtained directly from the spline interpretations of these 
numbers. We prove the main results in Section 4. 

2 definitions and notations 

Eulerian number, denoted here by A d ^ is the number of permutations in the 
symmetric group S d that have exactly k — 1 descents. An recurrence formula 
for A d ,k is 

A d:k +i = (k + + (d - k)A d ^ hk (2.1) 

with the boundary conditions 

4>,o = l,A d>0 =0,d> 0. 
An explicit formula for A d ^ is 

A d , fe = ^( d + 1 V-l) i (fc-i) d (2.2) 

which can be easily verified using the above recurrence. 

Refined Eulerian number A. d ,k,j is the number of permutations in the sym- 
metric group S d with k descents and ending with the element j . The recurrence 
formulas for A-d,k,j are [15] 

Ad+i,k,d-j+i — (k + l)A difeid _j + (d — k)A dl k-i,d-j, 

(2.3) 

A-d+i,k,d-j+i — kA dy k, d -j+i + {d— k + l)Ad t k-i,d-j+i- 
An explicit formula for A. d ,k,j is [15] 

Ad >fcJ = ^^(-^(fc-^-^k-i+ir 1 . 
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Descent polynomials, denoted by D d (t), arc defined as 



d 



k=0 



where D(d, n, fc) is the number of indexed permutations in S d with fc descents. 
The indexed permutation of length d and with indices in {0, 1, . . . , n — 1} is an 
ordinary permutation in the symmetric group Sd where each letter has been 
assigned an integer between and n — 1. Indexed permutations, or r-signed 
permutations, are a generalization of permutations (see jTTl [7] ) . We will follow 
the notation in [7J. The set of all such indexed permutations is denoted by S d l . 
The numbers D{d,n 1 k) satisfy a simple three-term recurrence [7], 



D(d,n,k) = (nk+l)D(d-l,n,k) + (n{d-k) + (n-l))D(d-l,n,k-l). (2.4) 
And also has an explicit formula |15j : 



We turn to the definitions of B-splines and Hcrmitc polynomials, then show the 
relations among Eulerian numbers, the refined Eulerian numbers, the descent 
polynomials and B-splines. 

Let p g R- 1 and L P (R) as usual denote the set 




k 



L p {R) = {/ : K -> C\f measurable, 




dt < oo}. 



For p = 2 and f,g e L 2 



(R) define the inner product 




and the norm 



11/11 := VUJ) 



making i 2 (R) to a Hilbert space. 
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For / <G L 1 (R) define the Fourier transform / A and the inverse Fourier 
transform f v as 



/+oo 
f(t)e- iut dt, 
-00 

n + OO 



oo 



2tt 



This definition can be extended to functions / € L 2 (K), see for example [5]. 
B-splines with order d, which is denoted as B d (-), is defined by the induction 

as 

Bt{x) 



1 if a: e [0,1), 
otherwise, 



and for d > 1 

B d = B x *B d _ u 
where * denotes the operation of convolution which is defined: 

/+oo 
f(t - V )g(y)dy 
-OO 

for / and g in L 2 (R). 

Then B d has the compact support [0, d] and is in C d_1 (R). A well known 
explicit formula for B d (-) is 

Bdix) = Jdhy. (9 (x - * } +" (2 ' 5) 

where the one-sided power function is defined by: 



x d x > 0, 

otherwise. 



And the recurrence relation is 



B d (x) = -^—B d - 1 (x) + ^-^B d - 1 {x-l). (2.6) 
d — 1 a — 1 

For all B d , d > 1, it holds S d G L X (R) f| L 2 (M) and 

^H = amc d (|) (2.7) 



5 



where sinc(t) denotes: 



{ 



sin t 



t = 0. 



sinc(t) : 



1 



(2.8) 



B-splines play important roles in approximation, computer-aid design(CAD), 
signal processing and discrete geometry. They have been well developed during 
the past few decades. For extensive monographs see [3J [H] . 

Hermite polynomials H n (x) with degree d, are set of orthogonal polynomi- 
als over the domain (— oo, +oo) with weighting function e~ x . The Hermite 
polynomials are defined by 



Being the limiting case of several families of classical orthogonal polynomials, 
they are of fundamental importance in asymptotic analysis. 

In this paper, we shall be concerned primarily with the development of var- 
ious asymptotic properties of Eulerian numbers, refined Eulerian numbers and 
the coefficients of descent polynomials, which are obtained directly from the 
spline interpretations of these numbers |15) . The theory of B-splines is a well 
developed area of applied numerical analysis and interpolation theory. Having 
related these numbers directly to B-splines [T5] , we can take advantage of many 
powerful spline techniques to derive various results of these numbers. It is pre- 
cisely this fact which motivates the need for more purely spline interpretations 
of these numbers to which we now turn. To state conveniently, we use [A J ]/(A) 
to denote the coefficient of X> in /(A) for any given power series /(A). 

Lemma 2.1. \W$ 




(2.9) 



and satisfy 




(i) 



A d , k 



d\-B d+1 (k); 
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(ii) D{d,n 7 k) = d\-n d -B d+1 [k + - ) ; 



(iii) A rf+1 , M _,- +1 = d\ ■ [X j ] UX + l) d B d+1 [k + A>0. 

3 Main results 

After the necessary definitions have been provided, we can come to the main 
results of this article. 

Theorem 3.1. Let be k € N then for d > k + 2 the sequence of the k-th 
derivatives of the B-spline converges to the k-th derivative of the Gaussian 
function, 



and 



d-voo 1 \12y d IV 12 2/J V2vr V 2 

where the limit may be taken point-wise or in L p (R),p 6 [2,oo). 

Sommcrfeld in 1904 [1 show that the Gaussian function can be approximated 
by B splines point-wise. In 1992, Unser and colleagues Q3] proved that the 
sequence of normalized and scaled B-splines B d tends to the Gaussian function 
as the order d increases in L P (R). A result due to Ralph Brinks [TB] generalize 
Unser's result to the derivatives of the B-splines. In this paper, we reprove the 
theorem and give the order of the convergence. 

Theorem 3.2. For x d = J^-x + we have 



^ = V^^ cxp (-Y) +0 (^ f ) (3 ' 3) 

L. Carlitz, D. C. Kurtz, R. Scoville and O. P. Stackelberg in [5] showed 
(|3.3|) with the help of the central limit theorem of probability theory that the 
expression on the right side of Q3.3P has the order 0(d~^~). Nevertheless, using 



the spline interpretation of Eulerian numbers, we obtain the same asymptotic 
forms of A n ^k with a much better convergence order 0(d~%). 

Theorem 3.3. For Xd = y^r^ + ^r-; we have 



^ DW „ J%)) = ^_^_ exp ^_(i±iLj +0 ( d -i) (3.4, 

Using the spline theory, we can also get the asymptotic representations of 
refined Eulerian numbers in terms of the Hermite polynomials. 

Theorem 3.4. Let x d = J~&$ (x - 1) + ^±1 then 

A d+1 , lxd] , d - j+1 = dl^ exp (-£-) g (±ti) "fli(*) + 



where H n (x) are the Hermite polynomials are defined by 

H n ( a: ) = (-l)» e Tr±_ e T-. (3.5) 

4 Proofs of the main results 



The proof of Theorem l3.1l is based on the following lemmata given in |T6! which 
shows an upper bound for standardized Sine functions. 



Lemma 4.1. 176] / For k, d € N and d > k + 2, iftere is a constant Cfc € K + suc/i 
i/iai /or 

GfcO) := Xr\[-i : i] (^)^^ + 7r fe |a;| fc exp(-x 2 ) 
where Xa{x) is characteristic function on set A, it holds 

G fc eL*(R).Vpe[l,oo) J 

and 



i ™ i 
7r a; 



7TX 



< G k {x). 



Proof. We start with the inequality 

Vx € [0, 1], sinc(x) < 1 — x 2 , 
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where the right-hand side term is the parabola that goes through the two ex- 
treme values of sinc{x) within this interval. We will first show that 

/ 2 \ d 

Vx G [0, Vd],sinc n (x/Vdj < 1-^- < exp (-a; 2 ) . (4.1) 
For this purpose, we define the positive function 

p(x) = exp (x 2 ) - . 

The derivative of p{x) is given by 

dp{x) _ -2exp(x 2 )x 3 (l-x 2 /d) d 
dx ~ d(l -x 2 /d) 

and is always negative for x £ [0, yd\. Therefore, the maximum of p(x) within 
the interval occurs at x = 



sup p(x)=p(0) = l, 



which proves (|4.ip . 

Second, for x > Vd, set 



Cfc := max 



fc + 2 

a 2 



We note that ir 



fcuiAs 



sinc(nx/ V") 



is also bounded by 



ir k \x\ k 



s/iic [ wx/Vd) < -^To- (4.2) 



The existence of Ck follows from the convergence of the sequence 

I n d-k+2 I 

\ / d>k+2 

Then one follows 

/ _\ d / _\ k+2 

v d \ ( v d \ 1 Cfc 



< — ~< 



We get a bound that is independent of d by noticing that 



Finally, we define the function Gk{x), which is independent of d, by suit- 
ably combining the right-hand sides of (14.11) and (14.21) . The function 7r fc |a;| fe • 
smc ^^2=1 is uniformly bounded from above by an L p {— oo,+oo) function 
Gfc(x) where p e [l,+oo). □ 



We have all the ingredients to prove our results now. 



Proof of Theorem \3.1[ Set 



L n (x) := d\n (sine i^J-J 



(4.3) 



Due to the symmetry, we may assume x > 0. By Taylor's theorem, for any 
x G [0, 1] and d € N, it holds 

/ x /l2~\ 1 x 2 1 „ ( 1 



and 



ln(l + x) = a; + G>(a; 2 ). 



Then for any a; € [0, 1] and d S N, it holds 

I -M-* ta ( 1 -T3 +0 (?))--T +0 (3 



Combing ()4.3[) and (|4.4[) . we have 

■ d /l2 \ 
sine —\—r — ex P 
I 2 V d ) 



•t) ( 1+0 G 



(4.4) 



(4.5) 



Furthermore, it holds = smc d (|), and S d e C" 1 - 1 ^)- Th e later yields 



for k < d- 1: 



(fc) 



(w) = i*uj k sinc a ( — 



Consequently, one obtains 



(A:) 



12* + 2 



i k uj k sinc d 



LJ 12 
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for d — > oo 



where the limit in equation (|4.6p follows from Lemma 14.11 and is meant to be 
taken point- wise. Clearly, B d k \uj) is bounded by the L p function Gk{x) defined 
in Lemma |4. II and this bound is independent of d. The use of Lebesgue's domi- 
nated convergence Theorem provides the L p convergence for P G [1, +00). The 
L q convergence with q 6 [2, +00] in the time domain follows as a consequence of 
Titchmarsh inequality, which states that for 1 < p < 2 and pT l + q^ 1 = 1, 
the Fourier transform is a bounded linear operator from L p {— 00, +00) into 
L q {— 00, +00). □ 

When k = 0, the Theorem 13.11 turns to 
Corollary 4.1. For d£N, the B-spline converges to the Gaussian junction, 



where the limit may be taken point-wise or in L p (R),p € [2, 00). 

Proof of Theorem \3.4\ We firstly prove Theorem 13.41 
Set 

p(A) = ^( J d )A d+1 ^^_ J+1 A-'" 

for A > — 1. Then p(X) is the polynomial of A with d degree. The lemm ^OI 
states that 

A d+lw+1 = d\ ■ [X j ] ((A + l) d B d+1 (k + )/(*), A >0. 
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which implies 

PW = E(') A "+i.M- J+ iA J = d\(X + l) d B d+1 (k + -i-) . (4.9) 

j=o ^ + s 

Taking the j-th derivative of both sides of equation (14. 9|) . it holds: 

,,/,,./.../■ = E(rf-j)!< 1 (fc + X ^ T )((A + l) d ) ( ^ ) ((A + ir 1 )('| A - i : 

■ =0 



J 



i=0 



(d-j+i)\ 



-j+i 

i=0 



Combining Theorem 13.11 with the definition of Hermit polynomials (|2.9I) , we 
obtain that 



□ 



Proofs of Theorem Iff. ,21 and Theorerr lS. 3\ By using Lemma 12.11 
(i) A d , k =dl-B d+1 (k), 



(ii) D(d, n, k) = d\ ■ n d ■ B d+1 [k+ — 

1 n 



and Corollary 14.1 



, d D d d\ 1 / x 2 \ (I 



consequently, for = i + we have 



= V W+T) cxp (-y) + 
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and 

□ 

Remark 4.1. We can also prove Theorem \3. 2\ by using Theorem \3.4\ By the 
combinatorial interpretations of Ad^k a>nd A.d,k,j> we have Ad+i.k = A^+i fc,d+i ■ 
Combining with Theorem \3.4\ one obtains 

^« = \/^TT5 <ix ' > (-?) +0 ^ § )' ( " 0) 
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